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Abstract. - Based on the Zimm-Bragg model we study cooperative helix-coil transition driven 
by a finite-speed change of temperature. There is an asymmetry between the coil^helix and 
helix^coil transition: the latter is displayed already for finite speeds, and takes shorter time than 
the former. This hysteresis effect has been observed experimentally, and it is explained here via 
quantifying system's stability in the vicinity of the critical temperature. A finite-speed cooling 
induces a non-equilibrium helical phase with the correlation length larger than in equilibrium. In 
this phase the characteristic length of the coiled domain and the non-equilibrium specific heat can 
display an anomalous response to temperature changes. Several pertinent experimental results on 
the kinetics helical biopolymers are discussed in detail. 



Biopolymers carry information which is embedded not 
only in the linear sequence of the monomers, but also in 
the conformational structures [1]. These structures are 
^ determined by the sequence, but they also adjust to en- 
i vironmental conditions. A pertinent example is the helix- 
' coil transition (HCT) , which denotes the disruption of the 
. ordered conformation — e.g., the a-helix of proteins or the 
' triple-helix of collagen — to form a disordered coil [1]. This 
order-disorder transition occurs when the temperature is 
raised or a chemical denaturant is added. 

The importance of helices in biopolymers motivated 
, many studies on thermodynamic and kinetic aspects of 
■HCT [1-6,9,10]. The basic model in this field was pro- 
\ posed by Zimm and Bragg (ZB) and successfully applied 
for describing the HCT both in [1-5] and out [6-12] of 
equilibrium. The virtue of the model is the simplicity of 
its ingredients: cooperativity and the free-energy prefer- 
ence to form a helix. Many studies devoted to the ki- 
netics of the HCT concentrate on the relaxation from a 
non-equilibrium state, which is prepared experimentally 
via, e.g., laser-induced temperature jump [10]. This setup 
is adequate for small globular proteins whose relaxation 
time is short. The situation is different for HCT in col- 
lagen [13-16], which has very long relaxation times to 
equihbrium (hours and days), so that the HCT has to 
be probed by necessarily finite-rate temperature changes, 
and the very description of the HCT has to be essentially 
kinetic [16] (the equilibrium HCT still exists at unrealis- 



tically long times [15]). Kinetic effects are also encoun- 
tered in HCT for DNA [17, 18] and for crystals of globular 
proteins, where each site of the crystal contains one pro- 
tein [20] . This ordered protein ensemble amplifies memory 
and hysteresis effects and allows their clear experimental 
identification [20]. 

These experimental set-ups call for a unifying theoreti- 
cal approach. Here we study the hclix-coil kinetics of ZB 
model driven by finite-speed temperature changes. We re- 
produce and explain several basic experimental findings, 
and also predict new effects. Conceptually, there is an 
even deeper interest in the kinetics of the basic ZB model, 
since the stability of many proteins does have both kinetic 
and thermodynamic aspects [21]. 

In the ZB model one assigns the spin variable Si = 1 
{si — —1) for the i'th helix (coiled) region of the polymer, 
and assumes the following free energy for the spins [1,2] 

F[S] = , - h{T)J2. S,, (1) 

^ — ^t— 1 ^ — ^i—l 

where s = (si,..sn), N is the total number of regions, 
and J > stands for the cooperative interaction [19]. The 
mechanisms of cooperativity for the main biopolymers (a- 
helices, DNA, collagen, etc) are reviewed in [1,2,14]. 

The term h{T) in (1) is the free energy difference of 
fast atomic variables in region z, calculated for a fixed 
value of the slow spin Si [27]. We thus assume time-scale 
separation [9,10]: the joint probability P(s,a) of fast (a) 
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and slow (s) variables factorizes as P(s)Peq(a|s), where 
the conditional probability Peq(a|s) is always (also for ki- 
netic processes) at equilibrium with the bath temperature 
T [27]. The dynamics of the spins is then governed by the 
free energy (1), where h{T) favors helix (coil) formation at 
low (high) T: h{Tc) = at the equilibrium HCT temper- 
ature Tc, while h{T) > {h{T) < 0) for T <T^{T> T^). 
Experiments and ah initio calculations are consistent with 
a linear change of h{T) in the vicinity of (a > is a 
constant and ks = 1) [1-4]: 

h{T) = a{n - T). (2) 

The interaction strength in the model is characterized 

by (T = e"^'^/"^^ [2]. In the highly-cooperative regime 
(7 <C 1 the equilibrium helix-coil transition resembles a 
real phase-transition, which combines the features of first- 
order (jumping order parameter) and the second-order 
(large correlation length) phase transitions [1,2,31]. 

The ZB model was originally proposed for describing the 
a- helix-coil transition in polypeptides and proteins [1,2]. 
Later on the equilibrium ZB model was successfully ap- 
plied to the duplex-coil transition in DNA [31] and to some 
aspects of the triplex-coil transition in collagen [14]. In 
the latter two cases the ZB model is regarded as a skeletal 
model producing important qualitative conclusions. The 
purpose of this Letter is to imderstand the basic physics of 
hysteresis and memory effects during the helix-coil transi- 
tion in terms of the ZB model. 

Now the system described by (1) interacts with a bath 
at temperature T. We assume that the elementary bath- 
driven process amounts to local disruption (or creation) 
of a single helix: Sj —Sj (spin- flipping), and that the 
dynamics is given by the following master equation [22]: 

^(«'^) = E,-=i iPi^j^tM^l^j) - P{s,t)w{sj\s)] , (3) 

where Sj = (si, .., Sj_i, — Sj, s^+i, .., sat), P{s,t) is the 
time-dependent probability of s = {si,..,Sn), P = dtP, 
and where the first (second) term in the RHS of (3) de- 
scribes the in-flow (out-flow) of probability to the configu- 
ration s due to spin-flipping. w(s]sj) is the transition rate 
Sj — > s which is standardly taken in the Glauber form [22]: 
w(sj|s) = ^[1 — Sj tanh(/?/ij)], where F is the relaxation 
frequency, (3 = 1/T and = h+J{sj-i+Sj+i) is the local 
field acting on sj; see (1). Thus the spin-flip is probable if 
it decreases F. This ensures relaxation to the equilibrium 
Peq[s] oc e-'^'^M for a constant T [22]. 

Note that the spin-flipping can occur anywhere in the 
chain. Thus, the studied kinetics of the ZB model differs 
from the zipper kinetics [8,11,12,17], where the disruption 
of the helix can occur only at the end-points of the chain. 
The zipper kinetics is expected to be valid for relatively 
short chains, undergoing relaxation from the completely 
helical to the completely coiled chains [8]. Here we con- 
sider long chains. 

Let us introduce the following averaged quantities: 

m{t) = {Si)t, €{t) = {SiSi+i)t, €2{t) = {SiSi+2)t, (4) 
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Fig. 1: The order parameter m versus dimensionless time vt 

under cooling-reheating with temperature speed v and param- 
eters (10), except that a = 1.2 x 10^^. Cooling starts at vt = 
and changes to reheating at vt — 0.4. The equiUbrium Tc is 
crossed for vt = 0.2 and t = 0.6. Thick curve: equilibrium rrieg, 
which changes between rrieq = — 1 (coil) and irieq = 1 (helix). 
Thick dashed curve: v/T = 2 x 10~^. Normal curves (from 
bottom to top): «/F = 5 x lO"'', 2 x 10^^ 10"^ The dashed 
counterpart of each normal curve refers to cooling from vt = 
till vt = 0.4, and then holding T constant. 

where (. . .)t is the average over P{s,t), while is the 
fraction of helical regions. For N ^ 1 (long chain) the 
boundary effects are neglected, all spins are equivalent, 
and (3, 4) imply 

r-im=-(l-ai)m+^ + ^e2, (5) 

r~-^ e = -2e-|- ai -I- (oo -h a2)m + aie2, (6) 
ai = /t+ — «;_, ao,2 = ±tanh(/3/i) + + K-,{7) 

where k± = ^tanh(/3/i ± 2/3J). Eqs. (5, 6) are flrst two 
equations of the infinite hierarchy of moment equations. 
For h = m = this hierarchy is exactly solvable [22] . For 
h ^ 0, there is no exact solution, and one has to rely on 
approximations [6,9]. 

The spin-temperature anzatz amounts to assuming that 
the probability P{s,t) has a locally Gibbsian form with 
two time-dependent parameters /3i and /32 [23,24]: 

P(s,i) ocexp[/3i(i)jV SiSi+i+/32(t)/iV.Si]. (8) 

This amounts to expressing the term €i in (5, 6) via e and 
m by means of equilibrium formulas [24] 

efc = TO2 + (l-rn2)l-'=(e-m2)^ fc = 2,3... (9) 

Thus we assume that the higher-order moments eft>2 re- 
lax to the local equilibrium (9), before m(i) and e{€) relax 
to equilibrium [7,23,24]. Eqs. (5, 6) with £2 given by 
(9) is now a closed pair of equations that reproduces ex- 
actly the equilibrium limit. Eqs. (5-7, 9) are consistent 
with the above-mentioned exact solution for h = m, = 0, 
since they correctly predict the diverging relaxation time 
for (3 J ::J> 1, as well as the transition from the expo- 
nential to a power- law decay [24]. The reliability of the 
spin-temperature anzatz is confirmed by its applications 
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Fig. 2: Susceptibility x{tjt') versus vt = t under cooling- 
reheating with parameters (10), except that a — 1.2 x 10~^ 
and v/T = 10~^. The temperature perturbation eTc9[T — t' + 
|]&[r' + f — t] has magnitude e = 0.001 and is centered at 
T. Normal curve: r' = 0.2; thick-dashed: t' = 0.3; dashed: 
r' = 0.41; thick: r = 0.6. Insert: the temperature profile 
C{vt) = T{vt)/Tc perturbed at vt = 0.2. 



in NMR/ESR physics [23,24]. For the derivations of (8) 
via the projection-operator method see [25]. 

The equilibrium m is obtained from putting to zero 



the LHS of (5, 6) with 63 given by (9): 



m, 



sinh(/3/i)[sinh^(/3/i) + e 



-4/3,71-1/2 



eq 



In the cooperative 



;-4/3„7 

leg — —1 (coil) to m,eq ~ 1 (helix) [1]; see Fig. 1. 



regime a = e <^ 1, rrieq shows a sharp transition 



from m 

Note that for /i = (right at T = T^) (5, 7) predict 
that the relaxation time of m{t) is ^(i-ai) ~ "^^^^ 
agrees with the known result [6] and shows that the relax- 
ation time is larger for more cooperative transitions and 
for lower temperatures. These aspects were numerously 
confirmed in simulations [12]. 

Now assume that the hath temperature T{t)/Tc = Co — vt 
decreases with speed v from a higher temperature T^Cq > 
Tc (at which the system was equilibrated) to a lower tem- 
perature TcQi < Tc- Then T{t) increases back to Tf,(o 
with the same speed (reheating): T{t)/Tc = Ci + 'vt- For 
identification of asymmetries in the system response, the 
reheating temperature profile is taken to be the mirror 
reflection of the cooling profile. 

For concretcncss the numerical solutions of (5, 6) are 
displayed for the dimensionless parameters [see (2)]: 

= 1.2, Ci = 0.8, a = 0.75, a = 3.34 x 10"*, (10) 

and various values of the dimensionless cooling-reheating 

speed v/T. The values of a and a correspond to the helix- 
coil transition in poly-7-benzyl-glutamate [3]. They are 
typical for other cooperative helix-coil transitions [4]. 
The order parameter m defines the helicity fraction 
. Fig. 1 displays the non-equilibrium m versus the di- 
mensionless time vt, as obtained from solving numerically 
(5, 6, 9) with the time-dependent temperature T{t). Fig. 1 
shows that for a small (but finite) speed v/T the transition 
helix— >coil during the reheating is more visible and takes 
shorter time than the reverse transition coil— >helix during 
the cooling. The same conclusion (not displayed) holds 
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when doing heating and then recooling. The symmetry be- 
tween liclix^coil and coil^hclix is recovered in the equi- 
librium limit 2^ <S 1- The asymmetry also disappears in 
the weakly-cooperative case it ~ 1. This asymmetry is an 
example of hysteresis and it was observed experimentally 
for highly-cooperative HCT in collagen [13-16], crystalline 
proteins [20] and DNA [18]. 

To gain a deeper understanding of the hysteresis, let 
us study in more detail the system's memory. Compare 
the cooling-reheating behavior of m(t) with the situation, 
where the temperature decreases till the lowest point and 
is then held constant; see Fig. 1. This cooling- holding 
scenario produces curves which are almost identical to the 
cooling- reheating curves, except at the vicinity of the equi- 
librium reheating transition, i.e., m{t) does not react on 
the reheating before the sign of h changes. Let us quantify 
the memory via susceptibility 



X{t,t') = lime^o[m(i) - m{t)]/e. 



(11) 



Here m{t) is obtained under the same cooling-reheating 
temperature setup, but e-perturbed at t': 

f{t)/Tc = Co-vt + se[vt + S- vt']0[vt' + S-vt], (12) 

where 9[t\ is the step function, and where the perturbation 
duration 2d is small but finite; see Fig. 2. The perturba- 
tion is designed such that at equilibrium, where m{t) is a 
function of the time-dependent temperature T{t), x(f,t') 
is non-zero only for v\t — t'\ < 6. In the regime where 
the above asymmetry is present, there are basically three 
scenarios for the behavior of x(^!^'); ^cc Fig. 2. i) A per- 
turbation introduced during the cooling in the vicinity of 
Tc is memoryzed and am,plified. This memory need not be 
monotonic: it revives once T(t) crosses Tc during the re- 
heating; see Fig. 2. ii) The same perturbation introduced 
during the reheating in the vicinity of Tc creates a stronger 
immediate response, but a weaker memory, as compared 
to the previous case, in) Outside of the vicinity of Tc the 
response resembles that in equilibrium: x(t; t') is maximal 
for t t' and quickly decays for |t — 1'| > 0; sec Fig. 2. We 
thus see how the hysteresis emerges from the unstabilities 
at T{t) ~ Tc during the cooling and reheating. 

Specific heat is a well-known indicator of the helix-coil 
transitions, which are observed via calorimetric methods 
[1]. Recalling the discussion before (2) and definitions (7), 
one can see that the energy of the spins is Nu{t) [27], with 



u{t) = -Je{t) - d^[ph{T)] m{t). 



(13) 



The origin of this formula should be clear from (1, 4). 
Note that the term — 5/3[/3/i(T)] is similar to the known 
equilibrium formula E = —d/3[PF] relating energy E to 
the free energy F. The specific heat c(t) is the response 
of u{t) to the temperature change [T = dT/dt] 

c{t) = u{t) /f{t) = -[Je + aTc m] /f. (14) 

The equilibrium specific heat Cgg is always positive and 
shows two sharp and symmetric peaks at T = T^; see 
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Fig. 3: Specific heat c versus the dimensionless time vt under 

cooling-reheating with parameters (10). Normal curve: v/T = 



X 10" 



v/r = 10-^. Thick curve: equilibrium Ceq- 

Fig. 3. For the non-equilibrium specific heat c we ace again 
the asymmetry between helix— >coil and coil— >hehx transi- 
tions: the peak of c during cooling is either absent or less 
visible than the one during the reheating. Now c{t) can 
be negative, i.e., the internal energy can decrease upon 
reheating; see Fig. 3. This is partially related to the re- 
sponse of m{t): Fig. 1 shows that when coohng changes to 
reheating, m{t) continues to increase due to its memory. 
In contrast, nieq decreases under reheating. The negative 
part of c is most pronounced for a finite cooling-reheating 
speed v; see Fig. 3. Note that a negative specific heat 
is met in glasses, within a different scenario that is also 
related to large relaxation times [26] . 

The kinetic transition temperature T can be related to 
the peak of the specific heat [15,16]; see Fig. 3. For not 
very small w, T is approximately a linear function of x = 
In [jfr] ; thus T is not susceptible to moderate changes in 
the cooling-reheating speed v. The same scaling of the ki- 
netic transition temperature was seen experimentally for 
the helix-coil transition of collagen [15]. For the parame- 
ters of Fig. 3 we obtained ('7c,r — Tc)/Tc = ac,RX + 6c, r 
for X G [—2, —3.5]. Here the lower indices C and R refer to 
the cooling and reheating, respectively, while or = 0.09, 
ac = -0.12, 6r = 0.32, and 6c = -0.41. We see that 
Tr > > To. 

The correlation function g{k,t) = {siSi^k)t " 
{si)t{si+k)t describes the spatial structure of fluctuations. 
Eq. (9) implies g{k,t) = [1 - m^it)] e-'^/^^*), where 



m = In 



1 - m^t) 
e{t)-rn^{t) 



(15) 



is the correlation length, or the cooperative unit length, 
which plays an important role in describing the coopera- 
tivity of the polymer structure [1]. The equilibrium ^eg; 
obtained from (15) by m{t) m^q and e{t) — * e^q, dis- 
plays two sharp peaks at T = T^; see Fig. 4. It directly 
relates to the cooperativity parameter: ^69(^ = ^0) = 
[1]. is small both below and above T^, since there are 
no much equilibrium fluctuations there. 

The non-equilibrium behavior of ^(i) is different; for fi.- 



Fig. 4: Correlation length ^ versus vt under cooling-reheating 
with parameters (10). Normal curve: v/V — i x 10^"^. Dashed 



Dashed curve: v/r = i x 10 . Dashed-thick curve: curve: 



v/V 



X 10" 



Thick curve: equilibrium 



nite speeds ^(f) increases during cooling, and its maximum 

is reached for the lowest T; sec Fig. 4. ^(t) is maximal at 
a certain finite speed of the cooling-reheating. Around its 
maximum it is larger than £,eq{T = T^)- Thus the equi- 
hbrium relation between step-like change of the order pa- 
rameter m and the correlation length is broken in kinetics: 
now m{t) does not show transition during cooling, but £^{t) 
is large. The reason for a large ^(t) is that the spin-spin 
interaction energy e{t) is close to 1 for both T < and 
T > Tc, while m{t) is far from ±1; see (15) and Fig. 1. 
Thus there are many helical and coiled domains, whose to- 
tal spin should sum to zero implying long correlations. Yet 
another interpretation of a large ^{t) is that due to mem- 
ory the system docs not sec the sign change of h{T), but 
it sees the lowering of T[t), which naturally increases its 
correlation length ^{t). Thus the finite-rate cooling plays 
a selective role siipprcssing one mechanism and activating 
another. In the non-cooperative case a ~ 1, £^{t) follows 
to the shape of ^eq, and ^(i) < for all temperatures. 

Domain lengths. Generally, the state of the studied lin- 
ear polymer is inhomogeneous and consists of helical and 
coiled domains. To characterize the domain lengths, we 

consider the probability V±{k) = ^11^=1^^^^ °f 

ing a hehcal {V+) or coiled {P-) domain of length k. 
The transfer matrix treatment of (8) leads to V±{k) = 



7±e 



-k/ri± 



where 7± and r]± do not depend on k. Thus, 



e-i^± =p^(2)/p±(l), or 



ri±{t) = In 



2 ± 2m{t) 
1 + e(i) ± 2m{t) 



(16) 



where 77+ and ?]_ arc, respectively, the characteristics 
length of the helical and coiled domains. Naturally, the 
equilibrium r]- eq is large in the coiled phase and decreases 
under cooling becoming small in the helical phase: see Fig. 
5. Likewise, rj+eq is large (small) in the helical (coiled) 
phase. Note that r]+eq{Tc) = r]-eq{Tc). In contrast to 
the equilibrium rj-eq, wc sec in Fig. 5 that r]-(t) can in- 
crease in time, if the cooling is not very slow. It decreases 
once the reheating starts. Thus ry_ reacts stronger on the 
decrease of temperature, than on changing the sign of h. 
Since during the cooling the order parameter m increases, 
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Fig. 5: The coiled domain length rj- versus vt under cooling- 
reheating witii parameters (10). Gray curve: v/T = ^ x 10~^. 
Thick-ashed curve: v/F = 10""^. Normal curve: v/T = i x 
10"^ Thick curve: equilibrium /^—^eg. 

we see that the number of helical segments increases, but 
the typical coiled domain gets larger. Fig. 5 also shows 
that ?]_ can behave non-monotonically with time. For 
very small speeds, r]- reproduces (with a delay) the shape 
of r]-eq; see Fig. 5. The behavior of ri+{t) under cooling- 
reheating is less interesting: it follows, in a delayed and 
weakened form, the shape of ?7_)_eg; see Fig. 6. (The same 
conclusion holds for r]^ during the heating from T < 
and then recooling.) Note that for finite speeds we can 
have r]- > 77+ at low temperatures; see Figs. 5, 6. 

Relations with experiments. We now discuss several 
experimental results on helical biopolymcrs that demon- 
strate clear signs of irreversibility and hysteresis. 

Admittedly, many experiments in polypeptides and pro- 
teins do not show visible signs of memory and hystere- 
sis (unless caused by irreversible aggregation), mainly be- 
cause the experimental temperature changes are too slow 
compared to relevant relaxation times [21]. 

However, there are susceptible experiments on crystals 
of globular proteins, which do see memory effects dur- 
ing various conformational changes [20]. Here each site 
of the crystalline structure contains one protein. Memory 
effects are amplified by this ordered structure and are visi- 
ble in experiments [20]. Along these lines, we present here 
several new experimental results. Note that irreversibili- 
ties and hysteresis effects here arc not associated with the 
crystalline structure per se [20]; they are related to helix- 
coil features proteins, while the crystalline structure only 
serves for amplifying these effects. 

Fig. 7 shows the denaturation process for crystalline 
protein Alcohol Dehydrogenase. The process is monitored 
via the change of the Young's modulus with temperature. 
Recall that the Young's modulus is defined as the ratio of 
the applied stress [prcssiirc] over the induced strain. The 
Young's modulus serves as an indicator of structural tran- 
sitions [28,29], since in the denaturated state it is smaller 
than in the native state. The Young's modulus of an Alco- 
hol Dehydrogenase sample was measured via analyzing the 
electrically excited transverse resonance vibrations of the 
sample, which is cantilcvered from one edge (another edge 
is free) [28,29]. The denaturation temperature of Alcohol 
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Fig. 6: The helical domain length versus the dimensionless 
time vt under cooling-reheating with parameters (10). Nor- 
mal curve: v/V = \ x 10"^. Thick curve: equilibrium T)-\-^e-q- 
Insert: the same with v/V = 10~^. 

Dehydrogenase (« 45.5° C) is identified via the sudden 
jmnp of the Young modulus; see Fig. 7. This agrees with 
the denaturation temperature obtained via calorimetric 
methods [30]. 

Fig. 7 shows that once the heating is substituted by 
re-cooling in the vicinity of the critical temperature, the 
system follows a different path (hysteresis), although the 
heating-re-cooling speed was rather small. Moreover, even 
though the heating has been changed to re-cooling, the 
Young modulus keeps on decreasing till 40° C due to 
the memory on the previous heating stage; see Fig. 7. 
These effects agree qualitatively with the theoretical re- 
sults found above via ZB model. 

In contrast to polypeptides and proteins, there are 
established experimental results concerning the hystere- 
sis and memory effects for DNA [17, 18] and collagen 
[13-16]. In both these biopolymers the helical state (du- 
plex for DNA and triplex for collagen) is stabilized by 
inter- molecular (i.e., inter-strand) interactions. 

Fig. 8 displays the experimental denaturation of an 
amorphous DNA film. The same effects of hysteresis and 
memory arc present here. For other experimental indi- 
cations of memory and hysteresis effects during DNA de- 
naturations see Refs. [18]. Ref. [31] critically assesses the 
applicability of the Zimm-Bragg model to the helix-coil 
transition in DNA, and finds that many experimental as- 
pects of this complex phenomenon are adequately reflected 
in the equilibrium Zimm-Bragg model. 

Ref. [14] investigates the equilibrium Zimm-Bragg 
model in the context of the helix-coil transition in col- 
lagen III. For this biopolymcr the end-points of the three 
strands are held together by disulfide bonds, which pre- 
cludes mismatches during the renaturation and makes pos- 
sible the application of the Zimm-Bragg model. Indeed, 
it was found that although experimentally the helix-coil 
transition in collagen III (as well as in collagen I) is always 
kinetic — the proper equilibrium regime is hardly reached 
within the experimental observation time — some impor- 
tant aspects of the phenomenon can be described within 
the equilibrium Zimm-Bragg model in quantitative agree- 
ment with experiments [14]. We thus expect that the ki- 
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Fig. 7: The Young modulus versus temperature for a rhombo- 
hedric crystall (direction [z]) each site of which contains Alco- 
hol Dehydrogenase protein. Molecular weight of this protein is 
80 kDa. The relative humidity is 97 % at 25° C. The water 
content at this temperature is 0.37 g of water per 1 g of dry 
protein. The speed of heating/recooling is 0.1° C per minute. 



Fig. 8: The Young modulus versus temperature for an amor- 
phous DNA film. Each DNA macromolecules (taken from a 
sturgeon mail) weights 1 000 kDa. The relative humidity is 95 
% at 25° C. The water content at this temperature is 0.42 g of 
water per 1 g of dry DNA. The speed of heating/recooling is 
0.1° C per minute. 



netics of the Zimm-Bragg model can describe the qualita- 
tive aspects of memory and hysteresis. 

In sum, based on the Zimm-Bragg model we studied ki- 
netics of the helix-coil transition driven by a finite-speed 
temperature change. We reproduced well-known experi- 
mental results on the hysteresis during the kinetic transi- 
tion and explained it by quantifying the process memory. 
We also predicted new scenarios of kinetic helix-coil tran- 
sition related to i) negative non-equilibrium specific heat 
accompanying the hysteresis; ii) correlation length becom- 
ing larger than in equilibrium. 

There is an increasing evidence that the characteristics 
of many important biopolymers is controlled by both ki- 
netic and thermodynamic factors [21]. For instance, the 
helix-coil transitions in collagen are normally kinetic, be- 
cause the equilibrium is not reached within reasonable 
times. Moreover, the kinetic helix-coil transition temper- 
ature of collagen for various organisms is close to their 
physiological temperature [15], since this kinetic transi- 
tion plays a role in achieving the flexibility of the collagen 
fiber [15]. Since our results indicate on new scenarios of ki- 
netic helix-coil transitions in the basic Zimm-Bragg model, 
they can be relevant for understanding the interplay be- 
tween the kinetics and thermodynamics in biopolymers. 

A. A. thanks Y. Mamasakhlisov for discussions. The 
work was supported by Volkswagenstiftung, ANSEF and 
SCS of Armenia (grant 08-0166). 

REFERENCES 

[1] K. Sneppen and G. Zocchi, Physics in Molecular Biology 
(Cambridge University Press, Cambridge, 2005). 

[2] H. Qian and J.A. Schellman, J. Phys. Chem. 96, 3987 
(1992). A.J. Doig, Biophys. Chem. 101-102, 281 (2002). 

[3] B.H. Zimm et at, PNAS, 45, 1601 (1970). 

[4] M. Go et al, J. Chem. Phys. 52, 2060 (1970). 

[5] M. Takano et al, J. Chem. Phys. 116, 2219 (2002). 



[6] G. Schwartz, J. Mol. Biol. 11, 64 (1965). D. Poland and 

H. A. Scheraga, J. Chem. Phys. 45, 2071 (1966). 
[7] H. W. Huang, Phys. Rev. A, 8, 2553 (1973). 

[8] T.R. Chay and C.L. Stevens, Macromolecules, 8, 531 
(1975). 

[9] A. Baumgartner and K. Binder, J. Chem. Phys., 70, 429 
(1979). 

[10] M. Takano et al, J. Chem. Phys. 118, 10312 (2003). 
[11] D. Porschke et al, Biopolymers, 12, 1313 (1973). 
[12] J. Bokkyoo and D.L. Weaver, J. Chem. Phys., 112, 4394 
(2000). 

[13] J.M. Davis and H.P. Bachinger, J. Biol. Chem. 268, 25965 
(1993). 

[14] J. Engel and H.P. Bachinger, Matrix Biol. 19, 235 (2000). 

[15] E. Leikina et al, PNAS 99, 1314 (2002). 

[16] C.A. Miles, Biopolymers, 87, 51 (2007). 

[17] V.V. Anshelevich et al, Biopolymers, 23, 39 (1984). 

[18] K.E. van Holde, Physical Biochemistry (Prentice Hall, NJ, 

1985). P.R. Bergethon, The Physical Basis of Biochemistry 

(Springer- Verlag, NY, 1998). 
[19] Note that the ZB model has several different appearances. 

All of them are equivalent to (1) for a long chain. 
[20] S.G. Gevorkian et al, Eur. Biophys. J. 34, 539 (2005). 
[21] D. Baker and D. Agard, Biochemistry 33, 7505 (1994). 

I. M. Plaza del Pino et al. Proteins 40, 58 (2000). 
[22] R.J. Glauber, J. Math. Phys. 4, 294 (1963). 

[23] M. Goldman, Spin Temperature and NMR in Solids 

(Clarendon Press, Oxford, 1970). 
[24] G.O. Berim and A.R. Kessel, Physica A 101, 112 (1980). 
[25] J. Ran and B. MuUer, Phys. Rep. 272, 1 (1996). 
[26] J. Bisquet, Am. J. Phys., 73, 735 (2005). 
[27] A.E. AUahverdyan and Th.M. Nieuwenhuizen, Phys. Rev. 

E 62, 845 (2000). 
[28] S.G. Gevorkian and V.N. Morozov, Biofizika 28, 944 

(1983). 

[29] V.N. Morozov and S.G. Gevorkian, Biopolymers 24, 1785 
(1985). 

[30] P.L. Privalov, Adv. Prot. Chem. 35, 1 (1982). 
[31] A. A. Vedenov et al Sov. Phys. Usp. 14, 715 (1972). 
M. Ya. Azbel, Phys. Rev. A 20, 1671 (1979). 



p-6 



